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ESTIMATION OF THE FORCES AND MOMENTS ACTING 
ON INCLINED BODIES OF REVOLUTION 
OF HIGH FINENESS RATIO 
By H. Julian Allen ' 

SUMMARY 

This report contains a discussion of the aerodynamic forces anri 
moments on inclined bodies of revolution. It is known that a simple 
potential flow solution for such bodies does not give results in good 
agreement with experiment. An approximate theory to allow for the effects 
of viscosity for such bodies is developed and applied. It Is shown that 
a simple allowance for viscous effects yields results in reasonable agree- 
ment with experiment for bodies of high fineness ratio such as would be 
used on missiles and supersonic aircraft. The methods developed are 
applicable at both subsonic and supersonic speeds. Some discussion of 
the probable effects of Reynolds number and Mach number on the forces nnfl 
moments on inclined bodies of revolution is included. 


INTRODUCTION 

Knowledge of the forces and moments on bodies of revolution has long 
been of interest in aeronautical engineering. The original interest per- 
tained to the characteristics of airship hulls, and one of the first logi- 
cal attempts at -understanding the nature of the flow field of these rela- 
tively long closed bodies was made by Max Munk (reference l) . Muhk demon- 
strated that on such closed bodies at a constant angle of pitch In straight 
flight and in a nonviscous fluid there occurred elemental forces along th e 
hull resulting from changes in the downward momentum of the fluid. Over 
the forward portions of the hull shown in figure 1 the downward momentum 
of the fluid must increase proceeding downstream because the apparent mass 
of the component flow normal to the axis of revolution increases due to the 
enlarging cross sections of the hull. Over this portion of the hull the 
reaction is upwardly directed for positive angles of attack. For bodies 
with parallel midsection, representative of the older airships, no cross 
force exists over these elements of the hull since there occurs no c hang e 
in momentum of the fluid as the air progresses along these sections of 
constant area. At the stern the contracting cross sections require a 
removal of momentum from the air stream and hence downwardly directed 
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elemental forces exist along the hull for positive hull angles of attack. 

It is shown in Munk’s work that for bodies of high fineness ratio the 
potential cross force per unit length fp at any station along the hull 
is given by 

f-n = (k 2 — k x ) q sin 2a (l) 

* dx — 

where S is the cross-sectional area of the hull; x is the distance along 
the hull from the bow, a is the angle of attack, 1 and k 2 and k x are, 
respectively, the transverse and longitudinal apparent mass coefficients 
for the body. The variation of k 2 — k x as a function of fineness ratio 
is given in figure 2. This cross force at small angles of-attack can be 
shown from the work of G. N. Ward, reference 2, to be directed ‘midway 
between the normal to the axis of revolution of the hull and the normal to 
the direction of motion of the hull (i.e., at an angle a/2). 

It is evident that for a closed body, such as an airship hull, at 
' a positive angle of attack the upwardly directed forces over the forward 
portion must be equal to the downwardly directed forces over the rear so 
that a pitching moment but no lift or drag results. In figures 3 and 4 
are shown a comparison of calculated and experimentally determined lift 
and pitching moments as a function of angle of attack for the hulls of 
the American airship ZR- 1 *- (U.S.S. Akron) (from reference 3) and the 
British airship R— 32 (from reference 4). It is seen that, contrary to 
the prediction of theory, a significant lift force exists at angle of 
attack. The pitching moment is reasonably well given by the theory, 
although the actual magnitude of the pitching moment is less, in general, 
than the predicted. 

Upson and Klikoff (reference 5) have compared the calculated and 
observed cross forces for several hull shapes, one of which is shown in 
figure 5< The discrepancy at the bow has been shown by Upson and Klikoff 
to be due to the bluffness of-the body. A discrepancy also exists between 
calculated and experimental characteristics over the remainder of the 
hull, the actual cross force over the rearward surfaces always being more 
positive for positive pitch than the potential theory predicts. It has 
long been recognized that this discrepancy results from the Influence of 
viscosity of the fluid. 

In recent years the determination of the cross forces on bodies of 
revolution haB again become of first Importance to the designers of super- 
sonic aircraft and missiles. These bodies in general differ from the usual 
airship hulls in two important aspects: First, the bodies are of higher 

fineness ratio, and, second, the bodies have a blunt stern or "base." 

■’’Throughout this report cross— force characteristics are considered in 
terms of the angle of pitch. It is clear that for a body of revolu- 
tion the arguments presented apply equally well to angles of yaw or 
any combinations of pitch and yaw. 
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It has been pointed out by Tsien (reference 6) that the older airship 
theories for the potential cross force are still applicable to such 
bodies even in supersonic flow for slender bodies at small angles of 
attack. For such unclosed bodies the airship theory would predict, in 
addition to a pitching moment, a lift and drag, but it is again to be 
expected that the failure to consider effects of viscosity will lead to 
important discrepancies between these calculated characteristics and experi- 
mental ones. 

It is the purpose of this paper to present an approx ima te analysis of 
the effects of viscosity on bodies of high fineness ratio. 


STMBOLS 

a speed of sound, feet per second 

A characteristic reference area of body for force and moment coeffi- 

cient definition, square feet 

Ap plan— form area, square feet 

c^ c section drag coefficient of a circular cylinder per unit length 
L in terms of its diameter 

Cp, body lift coefficient 

Cm body moment coefficient about an arbitrary axis at distance 
x m from the bow 

Cp body drag coefficient 

Cp a _Q body drag coefficient at zero angle of attack 

ACp increase in body drag coefficient above that at zero angle of 
attack 

D body drag, pounds 

Da=o body drag at zero angle of attack, pounds 

f cross force per unit length along the body, pounds per foot 

fp potential cross force per unit length along the body, po und s per 

foot 

f v viscous cross force per unit length along the body, pounds per 

foot 
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viscous cross force per unit length along the body of infinite 
length, pounds per foot 

longitudinal apparent mass coefficient 
transverse apparent mass coefficient 
actual body length, feet 

equivalent length of a blunt based body, feet 
lift, pounds 

moment about an arbitrary axis at distance x m from the bow, 
foot-pounds 

free-etream Mach number 

cross Mach number ( i . e component Mach number of the flow normal 
to the axis of revolution of the body) 

free-stream dynamic pressure, pounds per square foot 

radius of the body at any station x from the bow, feet 

maximum body radius, feet 

free— stream Reynolds number 

cross Reynolds number (i.e., Reynolds number based on the cross 
velocity) 

cross-sectional area of the body at any station x from the bow, 
square feet 

cross-sectional area of the base of the body, square feet 
free-stream velocity, feet per second 

cross velocity (i.e., component of the flow velocity normal to the 
axis of revolution), feet per second _ 

total volume of the body, cubic feet 

longitudinal distance from the bow, feet 

distance of the centroid of the plan-form area from the bow, feet 
station of the axis of moments, feet 

reference length used in definition of moment coefficient, feet' 


* r " - 
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a angle of attack, degrees or radians as ind icated 

T) ratio of the drag coefficient of a circular cylinder of finite 

length to that for a cylinder of infinite length 

v free— stream kinematic viscosity, feet squared per second 

p mass density, slugs per cubic foot 


THEORY 

In reference T, R- T. Jones, following earlier work of L. Prandtl in 
reference 8, considered the effects of viscosity on the flow over infin- 
itely long yawed cylinders and demonstrated that in the case of a laminsr 
flow the viscous effects may be considered by treating the flow across 
the cylinder axis independently of the flow along the cylinder. For cir- 
cular cylinders of infinite length the viscous force along the cylinder 
is simply that due to surface shear. The component flow across the cylin- 
der, however, introduces large cross forces due to separation of the flow. 
Jones has shown that the cross force on a yawed cylinder is accurately 
determined by considering the cross component of drag as may be seen in 
figure 6 taken from reference J. Although the demonstration of reference 
7 applies to laminar flows, it will be assumed to be applicable to turbu- 
lent flows as well. 


Consider, now, a body of revolution of high fineness ratio. It is 
again to be expected that the cross— force characteristics could be approxi— 
mately predicted by treating each circular cross section as an element of 
an infinitely long circular cylinder of the same cross-sectional area. 

With this assumption the local cross force per unit length due to viscosity 
fv^ would be given by 


fv 


CO 


2 ro i 0 



( 2 ) 


where r is the body radius at any station x from the bow, V c is the 
cross velocity, p is the mass density, and ca^ is the drag coefficient 
of a circular cylinder at the Reynolds number 


Rc 


2rVc 

V 


(3) 


and the Mach number 


M 


c 



w 


where, in addition, v is the kinematic viscosity, and a is the speed 
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of sound in the undisturbed stream. 

Since the cross velocity 

V c = V Q sin a (5) 

it follows that the viscous cross force becomes 

f v = Src^ q sin 2 a (6) 

where q is the dynamic pressure. The cross drag coefficient c^ c is 
that of a circular cylinder at the cross Reynolds number 


R c = 


2rVo 

V 


sin a = R q sin a 


( 7 ) 


and the cross Mach number 


M c = M 0 sin a (8) 

where Mq is the Mach number of the free stream. 

It is known (see appendix) that the drag coefficient of a circular 
cylinder of finite length is less than that for a cylinder of infinite 
length. A similar characteristic is to be expected as regards the viscous 
cross force for a body of finite length in oblique flow so that the vis- 
cous cross force will be less than that given by equation (6). It iB 
almost certain that the largest portion of the drag reduction due to 
finite length occurs at the ends of the cylinder. It will be assumed, 
however, that the reduction in drag for fineness ratio is the same for 
each element of a body of finite length so that in that case the viscous 
cross force becomes . 


f v = 2R rc dc <3 sin 2 a (9) 

where t} is the ratio of cross drag coefficient for the body of finite 
fineness ratio to that for a body of infinite fineness ratio. 

The integrated viscous cross force is then 

2r\q sins a J' rc^ dx 
o 

where l is the body length. 

In determination of the lift and drag characteristics at angle of 
attack, it should be noted that there also exists a viscous axial force 
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which is approximately the total drag at zero angle reduced hy the reduc- 
tion in axial dynamic pressure. The viscous axial force is then 

C-n q A cos 2 a (10) 

" 0=0 

where < -'D cl _q I s t ^ ra S coefficient at zero angle and A is the area 
upon which this coefficient is based. 

It is now assumed that the potential solution of Munk and the viscous 
solution may be combined to determine the cross— force distribution along 
the body and the integrated forces and moment on the body. The potential 
cross force per unit length acts at an angle a/2 from the normal to the 
free— stream direction, while the viscous force acts normal to the axis of 
revolution of the body so that the distribution of cross force in terms 
of the dynamic pressure is given, from equations (l) and (9), by 


f 

q 



a 

cos — + 
2 


q 


or 


L = (k-j-kiJ sin 2a cos | + 2 ^ 0 ^ sin 2 a 
The lift coefficient in terms of the reference area A 


(ID 


is given by 


_L_ 

qA 


where L is the total lift. Consideration of the potential cross force, 
the viscous cross force, and the viscous axial force then yields 


CL 


(k 2 -ii) sin 2a cos p- 
A 



2r\ sin 2 a cos a 

+ 

A 



dx — 


cos 2 a sin a 


( 12 ) 


The drag coefficient is obtained from the cross and axial forces as 


C D 


D_ 

qA 


where D is the total drag, or 
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C D = 


(k 2 -ki) sin 2a sin ^ AS 2t\ sin 3 a 


a" aa 

J dx 


dx + — / rc^ dx + Cr> 

J a c v a=0 


3 

cos a 


A 


( 13 ) 

The moment 1 coefficient about a given station ^ is dependent solely 
upon the cross— force distribution and is given by 


■'M 


C M “ 


(k 2 -k x ) sin 2a cos n 2 


AX 


I / as 
J dx 

o 


M 

qAX 

, 2t\ sin 2 a 

■*) ^ + jjT J rc d c ^ 


( 14 ) 


where X is a characteristic length for the evaluation of moment coeffi- 
cient. 


The method of the present report is clearly too approximate to justify 
the implied accuracy of the preceding equations.. It is considered that, 
in general, the following simplifications are warranted: 

1. Cosines of angles should be replaced by unity and sines of angles 
by the angles. 

2. The factor k 2 -k x should be replaced by unity. 

3. The viscous axial force (third) term of equation (12) may be 
neglected, while the corresponding term in equation ( 13 ) may be replaced 


'a=0 


Moreover the potential term integrals may be evaluated as 


and 



o 



— (z m -x) dx = vol - S b (l-x m ) 
dx 


where 


is the area of 


the body base and vol is the body volume. 


•The viscous cross— force term integrals may also be evaluated as 
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o 


and 




Ap ( Zjq— Z p ) 


where A^ is the plan— form area and Xp is the distance to its centroid 
from the how. y 


With the indicated changes equations (ll) to (lh) "become, respec- 
tively (with a in radians ) , 


and 


— = 2 ( a + 2 ti c j ra 2 

9 \dx/ c 

(11a) 


(12a) 

- (r)« 2 + v=4 c (x) a 3 

( 13 a) 


(I**) 


To determine the force and moment characteristics it is necessary to 
evaluate the coefficients ti and c& c . In the appendix, available data 
and some discussion of the coefficients are given. 


Comparison of Theory and Experiment 

Tests were recently completed at the Ames Aeronautical Laboratory 
of the high— fineness-ratio body shown in figure 7. The tests were con— • 
ducted in the subsonic speed range in the 12— foot pressure tunnel (refer- 
ence 9) , and in the supersonic speed range in the 6— by 6— foot wind tunnel 
(as yet unreported) . The test results afforded a good opportunity to 
compare the theory of this report with experiment. 

In none of the tests did the cross Mach number (given by equation 
(8) ) exceed 0.3 so that, aB may be seen from the appendix, no sensible 
compressibility effect exists, while the cross Reynolds number remained 
in the range for which c^ c is constant and equal to 1 . 2 . The coeffi- 
cient x\ was assumed to be dictated by the actual fineness ratio 
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_L_ = 3 • 9T2. a? 9.9 
2r 0 0.396 


for which 


t] =0.68 

With these "values and. the geometric parameters obtained from the 
shape equation (see fig. 7) the equations for lift coefficient, drag- 
coefficient increment, and moment coefficient my be determined from 
equations (12a) to (14a) as (with a in degrees) 


C T = i- = 0.019 a + 0.0025 a, 2 
L qA 

AC = = 0.00017 a 2 + 0.00004-3 a 3 

D • qA 

C M = — = 0.018 a + 0.00035 a 2 

M qAZ 

where A and l are the maximum cross-sectional area and the actual 
body length, respectively, and the pitching moment ie about a point 

x m = 0.704 Z 

which i6 the point about which the moments were considered in the experi- 
mental investigation of reference 9* The calculated characteristics of 
lift, drag increment, and pitching moment are compared with the experi- 
mental values in figure 8. It is seen that the theory of this report 
well predicts the lift and drag variation with angle of attack. The 
magni tilde of the experimental values of the pitching moment are lower , 
however, than the theory of this report would predict. 


DISCUSSION AND CONCLUDING REMARKS 

It is evident" from equations (12a) and (l4a) that the variation of 
lift and pitching moment of a body of revolution with angle of attack 
must, in the general case, be nonlinear since the potential cross forces 
due to the change of momentum of the fluid varies as the first power in 
a while the viscous cross force varies as the second power of the angle 
of attack. This is a particularly important aspect in the design of mis- 
siles since the guidance and control problem will be affected if nonlinear 
characteristics exist. It is of interest to note from equation (l4a) that 
the pitching-moment variation about the center of gravity for constant 
cross— force drag coefficient will be linear if the center of gravity is 
located at the centroid of the plan— form area. " 
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From the theory of this report, it appears possible to make the 
pitching-moment variation about the center of gravity of the body zero in 
one special case. This is the case of the body for which 


M = 2*r §£. r 
dx dx ~ 


that is, a conical body 


r fi x 

if the center of gravity is at the centroid of the cross area. 

Several other important expectations may be implied from the theory 
of this report: First, since the cross Reynolds number, which determines 

the cross 'drag coefficient, varies with angle of attack, it is possible 
for the cross drag coefficient to rather abruptly change (when R c varies 
with angle from 2 X 10 5 to 5 X 10 5 ) with angle of attack which might lead 
to rather erratic variations of the forces and moments with angle of attack. 
However, it is to be expected that the three-dimensional effects previously 
mentioned will reduce the indicated abrupt behavior. It iB also to be 
noted that wind-tunnel tests at lower scales would not necessarily show 
these peculiarities. Second, for bodies moving at high speeds the cross 
Mach number will increase from subsonic to Bupersonic values as the angle 
of attack is increased and the cross drag coefficient may vary in an 
erratic manner. This variation may again lead to corresponding variation 
in the forces and moment with angle of attack. Thus it iB clear that 
model tests of high-speed missiles should be performed over the whole 
speed range expected for the configuration if the model tests are to be 
Indicative of the true behavior. However, the fact that cross Reynolds 
number as indicated in the appendix is not important at Mach numbers above 
0.5 indicates that wind-tunnel tests on small-scale models at high super- 
sonic speeds should accurately predict the behavior of the full-scale 
configurations at and above the angle of attack for which the cross Mach 
number exceeds 0.5* 


Ames Aeronautical Laboratory, 

national Advisory Committee for Aeronautics, 
Moffett Field, Calif. 


APPENDIX 

The section drag coefficients of circular cylinders have been deter- 
mined for a wide range of Mach and Reynolds numbers by A number of experi- 
menters. A fairly comprehensive discussion of the drag phenomena of cir- 
cular cylinders is given in reference 10. 
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In figure 9 is shown the drag coefficient c* as a function of 

U-C . 

Mach number for circular cylinders of different sizes (corresponding to 
different Reynolds' numbers) . These values were obtained by W. F. Lindsey 
(reference 11)., John Stack (reference 12), T. E. Stanton (reference 13), 
and A. Busemann (reference l4) , as well as from some unpublished tests 
performed in the Ames 1— by 3— -l/2— : foot high-speed wind tunnel. 2 It will 
be seen that Reynolds number appears of significance only at low Mach 
numbers, so that for values of cross Mach number higher than 0.5 the 
curve of figure 9 may be expected to apply for all Reynolds numbers 
higher than about 10 2 . The variation of the drag coefficient c^ c with 
Reynolds number is shown in figure 10 along with some of the high sub- 
sonic drag characteristics shown in figure 9 and with the curves of E. F. 
Relf (reference 15) and C. Wieselsberger (reference 16) . Between figurtes 
9 and 10, the drag characteristics of circular cylinders as a function, of 
Reynolds and Mach number are fairly completely established. 

The position is not so fortunate with regard to t|, the ratio of the 
drag coefficient of circular cylinder of finite length to that of a cir- 
cular cylinder of infinite length, in that this ratio, to the author’s 
knowledge, has only been determined at one Reynolds number (88,000) and 
at a negligibly low Mach number (reference 10). These results are given- 
in figure 11 and correspond to the Reynolds number range for which 1.2 is 
the drag coefficient of the cylinder of infinite length. 

To obtain a rough estimate of the value of ^ at other Reynolds and 
Mach numbers, the following conjecture is given. The end— relieving effect 
for a cylinder of finite length must be primarily conveyed to other sec- 
tions through the low— velocity regions in the wake. Evidently the ratio 
of the spanwise length of the wake to the wake thickness would be the ratio 
which should determine q . The spanwise length of the wake will be approxi- 
mately the length of the cylinder, while the wake thickness will be nearly 
proportional to the product of the cylinder diameter and the drag coeffi- 
cient. It appears, then, that the value of q at Reynolds and Mach numbers 
for which ca c is not 1.2 would be the value of tj (from fig. 11) for an 
effective cylinder length— to-diameter ratio equal to the product of the 
actual length— to-diameter ratio and the ratio of the drag coefficient 1.2 
to the section drag coefficient at the Reynolds and Mach number in the case 
considered. 


g 

The 1— by 3— l/2-foot tunnel values (mainly useful in defining the trend 
at high subsonic values) were obtained using a rake of unshielded total- 
head tube and indicated drag coefficients about 15 percent higher than 
those obtained by others. It is believed that this effect was due to 
excessive angularity of the flow at the rake which would indicate incor- 
rectly high values. The values have been proportionately reduced to 
agree with Lindsey’s values and this proportionate reduction has been 
applied at all other Mach numbers . 
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It should, not he considered, that, because the section drag coeffi- 
cients have been given for supersonic Mach numbers, the equations devel- 
oped in the report are applicable at supersonic cross Mach numbers. The 
potential solution of Munk is derived on the assumption of an incompress- 
ible flow. This momentum solution, however, should be reasonably accurate 
up to cross Mach numbers of the order of 0.4. Lighthill (reference 17) 
has treated the problem of the inviscid cross force on bodies for the case 
in which the cross Mach number is not necessarily small. The solution 
obtained is given in increasing powers of the angle of attack. Lighthill*s 
solution, although numerically complex, may serve to replace the inviscid 
portions of the equations of this report. 
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Figure 2 inertia factors for bodies of revolution. 
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Figure 3- Lift and pitching- moment characteristics of a hut! model 
of the US.S. Akron. 
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Angle of attack , a, deg 

Figure 4. - Lift and pitching-moment characteristics of 
a hull mode! of the R-32. 
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Figure 5- Caicu/ated and experimental cross- force distribution on a mode! of the semi-rigid airship 
RS-t. 
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Angle between flow direction and the 
normal to the wire axis, deg 


Figure 6. - Variation of cross force on an oblique wire of 
circular cross-section . 
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Figure 7.- Body of revolution employed for the example of this report. 
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Lift coefficient 
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(a) (. iff characteristics. 


Figure 8 Calculated and experimental aerodynamic characteristics 
for the example body of revolution. 
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(b) Drag characteristics. 


Figure 8.- Continued. 


COHESEKEHI^ 



Pitching-moment coefficient 


24 



Angle of attack, a, deg 


(c) Pitching-moment characteristics. 


Figure 8- Concluded. 
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Figure 9. — Drag coefficients of circular cylinders of various sizes as a function of Mach number. vn 
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